Can hadronic rescattering explain the "jet quenching" at RHIC? 
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Recent RHIC data have shown novel nuclear modifications of moderate to high pr particle produc- 
tion in central Au+Au collisions, including a suppression of hadron production and a disappearance 
of back-to-back hadron pairs. In this paper, we investigate whether final-state hadronic interactions 
of the jet fragments can reproduce the RHIC data. We find that hadronic rescattering can account 
for the disappearance of back-to-back hadron pairs, but cannot reproduce other features of the 
RHIC data. 
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Recent data on the production of high px hadrons in 
central Au-|-Au collisions at RHIC indicate novel nuclear 
effects. Above pr — 5 — 6 GeV/c, the yield of hadrons is 
suppressed by a factor of w 5 compared to what would 
be expected from an incoherent superposition of inelastic 
nuclcon-nuclcon collision 0, Q . Data on two-particle az- 
imuthal correlations show similar near-angle jet-like cor- 
relations in central Au-|-Au and p+p collisions 0, 
The back-to-back dihadrons indicative of dijet produc- 
tion, however, are absent or greatly suppressed in the 
most central Au+Au collisions. In addition, the produc- 
tion of high pt hadrons shows a strong azimuthal cor- 
relation with respect the reaction plane ("elliptic flow") 
i 

Taken together, these experimental data are thought 
to result from a novel nuclear effect known as jet quench- 
ing. Recent measurements from d-|-Au collisions do not 
show the same behavior as central Au-|-Au collisions [3, 
so the modification of moderate to high px hadron pro- 
duction observed in central Au-|-Au collisions are thought 
to arise primarily from the interaction of fast partons or 
their fragmentation products with the dense medium pro- 
duced in central collisions of heavy nuclei. The goal of 
the current work is to investigate whether the data can 
be explained entirely in terms of the hadronic interac- 
tions of jet fragmentation products with a dense hadronic 
medium. 

Fast partons traversing a dense gluonic medium are ex- 
pected to lose energy and acquire transverse momentum 
relative to their original direction of propagation Ji,.i,^J. 
This energy loss is due to radiative induced gluon emis- 
sion. The rate of energy loss is proportional to the gluon 
density of the medium traversed. The previously men- 
tioned RHIC data have been described quite successfully 
by convoluting expected parton production rates, parton 
energy loss in an expanding dense medium, and parton 
fragmentation. 

To quantify the nuclear matter effects on high px par- 



ticle production the nuclear modification factor Raa is 
constructed, 

TAAd^a^^/dpTdTj ' ^ > 

where the nuclear overlap function TAA—{^hin) /cTi^ei 
from a Glauber calculation accounts for the nuclear col- 
lision geometry. In the absence of nuclear matter effects, 
Raa should approach unity at moderate pr (2-3 GeV/c). 
In central Au-f Au collisions at RHIC, Raa ~ 1/5. This 
behavior was predicted by models that incorporate par- 
tonic energy loss in a dense gluonic medium. A recent 
paper, however, pointed out that this suppression of sin- 
gle inclusive particle production at high px can also be 
explained qualitatively by assuming that partons frag- 
ment inside a dense hadronic medium The basic 
feature of the single inclusive data {Raa ~ 1/5) is re- 
produced by this model, although it may fail to describe 
the exact pt dependence of the suppression. Adopting 
the notion of a colorless pre-hadron allows for a better 
description of the px dependence of Raa 'i3| . With suit- 
able modifications to the initial state nuclear effects it is 
likely that this hadronic rescattering model could give a 
reasonably qualitative description of the single inclusive 
data, including the exact pr dependence of Raa- 

The ability to describe the single inclusive data both 
in terms of parton energy loss in a dense gluonic medium 
and hadronic rescattering and energy loss in a dense 
hadronic medium is not surprising. Any many-body cal- 
culation based on the Boltzmann equation yields iden- 
tical single particle distributions under the substitution 
a — > Aa and p p/A where a is the two-body scat- 
tering cross-section, p is the density of scattering cen- 
ters, and A is an arbitrary number. Thus, a single par- 
ticle observable is unable to distinguish between a dense 
medium with small scattering cross-sections and a more 
dilute medium with larger scattering cross-sections pdj . 
In general, QCD energy loss cannot be described using 
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the Boltzmann equation as QCD energy loss processes 
may be coherent. Nonetheless, the uncertainties in the 
dynamical evolution of a heavy-ion collision coupled with 
the uncertainties in nuclear effects in parton production 
from cold nuclei make it difficult to distinguish between 
the partonic and hadronic energy loss scenarios using sin- 
gle particle data alone. 

Fluctuation observables are not invariant, however, un- 
der the substitution cr — * Aa and p p/A. For A > 1, 
fluctuations will be increased. Fluctuations can be mea- 
sured via many observables, but are most easily quanti- 
fied in terms of two particle correlations. The production 
of jets in a heavy-ion collision can be thought of in terms 
of a local fluctuation in high pT particle production. In 
this paper, we investigate what a hadronic rescattering 
interpretation of the RHIC single inclusive particle pro- 
duction data would predict for two-particle azimuthal 
correlations. We find that the hadronic rescattering pic- 
ture is unable to simultaneously describe the RHIC data 
on inclusive particle suppression and high-p-r azimuthal 
correlations. The paper is organized as follows. We first 
compare RHIC directly to the PYTHIA event generator 
and find that this model reasonable describes the fea- 
tures of moderate to high px inclusive particle production 
and two-particle azimuthal correlations. We then intro- 
duce a hadronic rescattering model that will be used to 
model a dense hadronic system into which we will em- 
bed fragmentation products from PYTHIA events. We 
then merge the PYTHIA jet events with the rescatter- 
ing model, and study the propagation of PYTHIA frag- 
mentation products in our dense hadronic medium. We 
find quantitative disagreement with the RHIC data, par- 
ticularly due to the copious resonance production from 
hadronic rescattering that should be manifest in the two- 
particle azimuthal correlations but is not observed in the 
experimental data. 



I. THE PYTHIA EVENT GENERATOR: 
COMPARISON TO RHIC DATA 

In order to model hard scattering processes and frag- 
mentation, we use the PYTHIA event generator 
This model convolutes measured parton distribution 
functions, elementary parton-parton scattering cross- 
sections and a phenomenological model of jet fragmen- 
tation. In Figure ^ we compare the invariant pr spec- 
tra for pions produced running PYTHIA with its stan- 
dard settings with the recently measured PHENIX p + p 
data 13] . For our studies here we are only interested 
that the shape of the PYTHIA spectra is similar to the 
data, so the PYTHIA spectrum is normalized to match 
the real data at = 4 GeV/c. As seen, the PYTHIA 
event generator with default settings does a reasonable 
job in describing the measured RHIC data between 4 and 
8 GeV/c. 

Our main goal in this paper is to investigate azimuthal 
correlations among highp^ hadrons. In Figure|21we com- 
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FIG. 1: Comparison of the invariant spectra from RHIC p + p 
data [13] and PYTHIA. The PYTHIA spectrum is arbitrarily 
normalized. 



pare the azimuthal correlations measured by the STAR 
collaboration 4] to those produced by the PYTHIA event 
generator. The PYTHIA calculations use identical kine- 
matic cuts as the STAR collaboration (pseudorapidity 
rj < 0.7). Events with a high px trigger hadron with 
4 < PT < 6 GeV/c are found. We then calculate the az- 
imuthal separation of other hadrons in these events with 
PT > 2 GeV/c. The resulting azimuthal distribution is 
normalized to the number of trigger hadrons. In both 
the STAR data and the PYTHIA calculations there are 
strong azimuthal correlations near A(j> « and A(j) w tt. 
The correlated back-to-back hadrons arise from the frag- 
mentation of back-to-back dijets. 
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FIG. 2: Per trigger hadron relative azimuthal distributions 
for STAR data ^ compared to PYTHIA calculation with 
and without resonance decays. 

The correlated small-angle pairs can be produced via 
two mechanisms. During jet fragmentation, the pro- 
duced hadrons will be collimated thereby producing az- 
imuthally correlated particle pairs. Such correlated pairs, 
however, can also be the result of the decay of resonances. 
In Figure |21 we investigate the role of resonance decays 
in the PYTHIA calculations. Within this model, sta- 
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ble hadrons (tt, K^p) are produced directly during string 
fragmentation and also during the decay of unstable res- 
onances produced during string fragmentation. In Fig- 
ure |21 we show two different PYTHIA calculations. The 
standard PYTHIA setting includes decay of unstable res- 
onances and we see a strong near-angle and back-to-back 
correlation peak. We also ran PYTHIA in a mode where 
no resonances produced during fragmentation were al- 
lowed to decay. In this mode near-angle and back-to- 
back correlations are also observed. These correlations, 
however, are greatly suppressed. It is important to note, 
however, that the suppression factor is identical for the 
near-angle and back-to-back hadron pairs. Since the 
back-to-back hadron pairs cannot arise from low mass 
resonance decay (the invariant mass of the parent would 
need to be greater than 5 GeV) we conclude that the 
near-angle correlations observed in PYTHIA arise pri- 
marily from the jet fragmentation and are not due to 
resonance decay. Turning on the resonance decays leads 
to the observation of more charged hadrons per event 
and hence increases the strength of the near-angle and 
back-to-back azimuthal correlations. Analogously, if we 
were to include tt*^ in the construction of the azimuthal 
distribution, we would see a trivial 50% increase in the 
near-angle and back-to-back correlation strength. The 
increase of the azimuthal distributions is due primarily 
to the observation of the extra hadrons from resonance 
decay and not from the correlations induced by these de- 
cays. 

It should be noted that the absolute per-trigger yield 
of associated particles differs between PYTHIA and the 
STAR data. In the case of PYTHIA with decays, the per- 
trigger yield of associated particles is larger for PYTHIA 
compared to the STAR data. In the case of PYTHIA 
without decays, the per-trigger yield of associated par- 
ticles is smaller than the STAR data. No attempt was 
made to adjust the PYTHIA fragmentation settings to 
rectify this difference. These azimuthal distributions are 
normalized to the number of "trigger hadrons" , and these 
mostly arise from the case where the leading hadron car- 
ries a large fraction of the parton momentum. This is the 
region where the experimental data on parton fragmenta- 
tion functions is most uncertain, and PYTHIA fragmen- 
tation functions are tuned to experimental data. Thus 
any discrepancy between measured dihadron correlation 
data and the PYTHIA calculations is not of fundamental 
concern. 



II. A SIMPLE HADRONIC RESCATTERING 
MODEL 

The hadronic rescattering model JJjJ that we will use 
to understand high pT hadron propagation in a dense 
hadronic medium has been described elsewhere and was 
shown to reproduce many of the low pr experimental ob- 
servables at RHIC including transverse momentum dis- 
tributions for TT,K,p, U2(pt), and the Hanbury-Brown 



Twiss radii Il5|. We briefly review some of the features of 
the model below. 

Rescattering is simulated with a semi-classical Monte 
Carlo calculation which assumes strong binary collisions 
between isospin-averaged hadrons. Relativistic kinemat- 
ics is used throughout. All calculations are made to sim- 
ulate RHIC-energy Au-|-Au collisions in order to compare 
with the results of RHIC data. 

The initial stage of the rescattering calculation em- 
ploys simple parameterizations to describe the initial mo- 
menta and space-time of the hadrons. The initial mo- 
menta are assumed to follow a thermal-like transverse 
momentum distribution for all particles, 

{l/mT)dN/dmT = CTOT/[exp (mr/T) ± 1] (2) 

where T is a "temperature parameter" , and a Gaussian 
rapidity distribution for mesons, 

dN/dy^Dexp[-{y-yof/{2ay'')] (3) 

where ay is the rapidity width. Two rapidity distribu- 
tions for baryons have been tried: 1) flat and then falling 
off near beam rapidity and 2) peaked at central rapidity 
and falling off until beam rapidity. Both baryon distribu- 
tions give about the same results. The initial longitudinal 
particle hadronization position {zhad) and time {thad) are 
determined by the relativistic equations, 

Zhad = Thad sinh y; thad = Thad cosh y (4) 

where Thad is the hadronization proper time. From Equa- 
tions 2 and 3, it is seen that longitudinal invariance is 
not assumed in the initial conditions for the present cal- 
culations. Calculations were carried out using isospin- 
summed events containing at freeze-out about 5000 pi- 
ous, 500 kaons, and 650 nucleons (A's were decayed). 
The hadronization model parameters which were found 
to reproduce the RHIC data were T = 300 MeV, (7^=2.4, 
and Thad=^ fm/c. 

We now compare RHIC high px data to calculations 
from the pure rescattering model. In other words, for 
these calculations hard scattering and fragmentation is 
not implemented. In Figure 13 we compare the PHENIX 
p + p tt" + AT data to the pure rescattering model 
calculations for Au-|-Au collisions with impact parame- 
ter 6 = 4 fm. For high pT 7r° production at RHIC, the 
spectral shape is similar in central Au-|-Au and p-l-p col- 
lisions. We are only interested in spectral shapes, so the 
rescattering model calculations are arbitrarily normalized 
to the PHENIX data at pt = 4 GeV/c (The integral of 
the pt spectra and hence the total yield is used to adjust 
the initial conditions of the calculation). We see that the 
rescattering model produces copious high px hadrons and 
exhibits a power law behavior. The rescattering model 
over-predicts the yield of high px hadrons. This is prob- 
ably due to the treatment of the densest stage of the 
collision in terms of 2 ^ 2 binary hadronic collisions 
whereas a more correct treatment would model this stage 
using many-body collision dynamics. It has been verified 
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numerically that this approximation does not affect the 
bulk dynamics of the collision 11 j as measured by low 
Pt observables. 
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FIG. 3: Comparison of the invariant spectra from RHIC p + p 
data and a pure rescattering model calculation for Au+Au 
collisions with b=4 fm. The rescattering model spectra is 
arbitrarily normalized. 

The pure rescattering model calculations also show 
near-angle (but not back-to-back) azimuthal correlations. 
Figure ^compares the measured azimuthal distributions 
for p-t-p collisions Q and the azimuthal correlations pro- 
duced in the rescattering model calculation for Au + Au 
at 6 = 4 fm. The near-angle correlations in the rescatter- 
ing model come from the decay of resonances. We have 
confirmed this by considering rescattering particles that 
do not come directly from a resonance decay, in which 
case no azimuthal correlations are observed. These cor- 
relations are both stronger and broader than those mea- 
sured in the p-l-p data. We have fit the rescattering model 
near-angle azimuthal distribution with a Gaussian and 
get a width of = 0.39 ± 0.04. In contrast STAR has 
measured the widths of the near-angle azimuthal corre- 
lation in p-l-p, d-f Au, and central Au+Au collisions and 
gets a width in all cases of ct w 0.2 for these kinematic 
cuts. The large discrepancy in the widths of these corre- 
lations indicates a fundamentally different origin. In the 
rescattering model, the correlations come from resonance 
decay, whereas for the RHIC data the correlations seem 
to arise primarily from jet fragmentation. 

The near-angle azimuthal correlation strength is re- 
lated to the yield of resonances at pt ~ Pt (trigger) + 
Pt {associated). The fact that the pure rescattering 
model produces an invariant px spectrum of hadrons and 
resonances that is harder than that measured in the data 
leads to a somewhat artificial increase in the near-angle 
correlation strength. This model is certainly incorrect 
for calculating moderate to high pT particle production 
in Au+Au collision at RHIC, so we chose to make no at- 
tempt to correct the model for its possible shortcomings. 
The rescattering model will only be used as a rough esti- 
mate of the environment that a jet fragmentation prod- 
uct would see if it was able to fragment inside such a 



dense hadronic medium. Nonetheless, we note that this 
pure rescattering model does produce azimuthal corre- 
lated particle pairs due to resonance decay, and that these 
correlations are broader and stronger than those seen in 
the real data. 
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FIG. 4: Per trigger hadron relative azimuthal distributions 
for STAR p + p data compared to a pure rescattering model 
calculation for Au + Au at fe = 4 fm. The blue line shows a 
Gaussian fit to the azimuthal distribution from the rescatter- 
ing model. 



III. THE PROPAGATION OF PYTHIA TEST 
PARTICLES IN THE RESCATTERING MODEL 

PYTHIA events are embedded into these rescattering 
events, with the production vertices chosen to account 
for the initial nuclear overlap geometry. The embedded 
PYTHIA events are required to have a high pT hadron 
with pt> 3 GeV/c. The particles coming from these 
PYTHIA events are tagged as such, and any resonance 
produced via the interaction of a PYTHIA particle is 
tagged as a "PYTHIA resonance". The decay products 
of these PYTHIA resonances are tagged as PYTHIA par- 
ticles, and in the plots that follow we only look at the 
final-state PYTHIA particles. Thus, the artificially large 
high Pt particle production seen in the pure rescattering 
model does not invalidate our subsequent calculations. 

The main theoretical uncertainty in these studies is the 
space-time development of hadrons in parton fragmenta- 
tion. In the standard PYTHIA fragmentation scheme, 
no attempt is made to describe the space time aspects of 
jet fragmentation. There are many theoretical models of 
how jet fragmentation develops, but little experimental 
support for any single idea. In the conventional QCD 
model, hadrons from parton fragmentation are formed 
at time t = ER^ 16], where E is the parton energy and 
R is hadronic size. Thus, very high pT hadrons would 
necessarily form well outside the medium produced in 
a Au+Au collision at RHIC. In another picture color- 
less pre-hadrons are formed very early ^3 in which case 
hadrons would form deep inside the medium produced 
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at RHIC. The pre-hadrons form instantly in the limit 
where — *■ 1, where Zh is the fraction of the parton 
momentum carried by the leading hadron. One argu- 
ment against such a hadron formation picture at RHIC 
is that hadrons should not exist in a quark-gluon plasma. 
The instantaneous formation of a colorless pre-hadron is 
based on the notion of vacuum energy loss and motivated 
by the string picture. In a quark-gluon plasma, Debye 
screening leads to the modification of the string tension 
and would reduce the rate of vacuum energy loss. 

Here we adopt a naive picture of hadron formation, 
investigate its consequences, and rule it out along with 
any similar fragmentation picture. In these calculations, 
the hadrons from the PYTHIA events initially form at 
proper time, r, and at (z,t) as given by Equation 3 with 
Thad = T = 1 fm/c, to be consistent with the rest of 
the calculation. In the transverse direction, the (x,y) 
position of the embedded PYTHIA events are randomly 
selected using the nuclear overlap model. The individual 
hadrons in each PYTHIA event are further smeared in 
a circle of radius 1.5 fm. This radius is determined by 
the condition that it is large enough such that hadronic 
rescattering within a PYTHIA event is small. 

Figure[Slshows the nuclear modification factor Raa for 
pions from these calculations as well as PHENIX Au+Au 
data Ij at a similar impact parameter. For the model 
calculations, Raa is defined as: 



(PYTHIA tagged from rescattering events) 

^^^(^^^ ^ A( PYTHIA p -hp) ■ 

(5) 

With this definition, we consider only the modification of 
the spectra of the embedded PYTHIA events, and do not 
consider the contribution to particle production coming 
from the rescattering model. Because we inject triggered 
PYTHIA events into the rescattering model, we only con- 
sider Raa for pt > ^ GeV/c. 

Figure |31 shows similar features to the calculations of 
Gallmeistcr ct al. 9] . The absolute value of Raa is well re- 
produced in the moderate region {j)t ~ 4 GeV/c), but 
Raa trends upwards with increasing , a trend not seen 
in the experimental data. The pr dependence arises from 
our treatment of the space-time development of hadrons 
in jet fragmentation. Higher pT hadrons are not able to 
interact until proper time r and they thus travel a dis- 
tance of 7T in the laboratory frame (here 7 is the Lorentz 
factor). 

Figure|Blshows the two-particle azimuthal distributions 
from the STAR collaboration Q and from the rescatter- 
ing events with embedded PYTHIA jets. For the model 
calculations, we only look at hadrons produced directly 
from the PYTHIA event or produced via the decay of a 
PYTHIA particle induced resonances. In the model cal- 
culations all hadrons are considered (charged and neu- 
tral) whereas for the STAR measurements only charged 
hadrons are used. This should lead to « 50% stronger 
correlations in the model compared to the real data. 

As pointed out by STAR, the correlations near A(j> « 
are indicative of jet fragmentation, and are very simi- 
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FIG. 5: Nuclear Modification Factor Raa for pions. The open 
squares show PHENIX Au + Au data for the 10% most central 
collisions, and the solid circles show the results of the rescat- 
tering model with embedded PYTHIA events. Only particles 
tagged as originating from PYTHIA or a PYTHIA-particle 
induced resonance are included in the rescattering-|-PYTHIA 
calculation. 
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FIG. 6: Per trigger hadron relative azimuthal distributions for 
STAR data compared to PYTHIA+rescattering model calcu- 
lations. Only particles tagged as originating from PYTHIA 
or a PYTHIA-particle induced resonance are include in the 
rescattering-|-PYTHIA calculation. 



lar in magnitude and width for proton-proton and cen- 
tral Au+Au collisions. The azimuthally back-to-back di- 
hadron pairs (A^ « tt), indicative of dijet production, 
are present in the proton-proton collisions but absent in 
the most central Au-f An collisions. This feature is well 
reproduced in the rescattering model -I- PYTHIA jet cal- 
culations. 

Focusing, however, on the region A0 « reveals a 
striking difference between the rescattering-|-PYTHIA 
calculations and the STAR data. While the STAR 
data shows little or no modification of the correlation 
structure in central Au+Au collisions compared to the 
proton-proton reference, the rescattering model leads to 
a broadening and enhancement of the near-angle corre- 
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lations. This is exactly the structure observed in the 
pure rescattering model calculations (where no PYTHIA 
jets were embedded into the events). These correla- 
tions are due primarily to the resonances produced dur- 
ing the rescattering stages of the collisions, and do not 
come from the initial jet fragmentation. We thus con- 
clude that the near-angle correlation structure observed 
by STAR in central Au-|-Au collision is due to jets that 
fragment outside the medium. If the jets were to frag- 
ment inside a dense hadronic medium they would pro- 
duce resonances that lead to broad azimuthal correla- 
tions. The broad resonance induced correlations seen in 
the PYTHIA-|-rescattering calculations are not seen in 
the central STAR Au-|-Au data. 

IV. CONCLUSIONS 

The failure of this model to describe in detail all as- 
pects of high pt particle production at RHIC is not sur- 
prising. It is unlikely that a 2 — > 2 hadronic scattering 
description is valid in the early stages of a RHIC colli- 
sion. For low Pt observables such as elliptic flow and 
HBT this model does seem to work. At high pT, how- 
ever, we show that this model does not work and we 



feel that no purely hadronic model can describe the "jet 
quenching" observables seen in central RHIC collisions. 
The nuclear suppression factor will have an unwanted pr 
dependence. In addition, the unavoidable copious reso- 
nance production will lead substantial modification in the 
two-particle high pT azimuthal correlations. In particu- 
lar, we find that resonance decays lead to broader high 
Pt azimuthal correlations, and there is no evidence in the 
RHIC data for large broadening of the azimuthal corre- 
lations. In general, we find that only correlation and 
fluctuation measurements, of which high pT azimuthal 
correlations are one example, are able to distinguish be- 
tween a purely hadronic description of RHIC data and an 
interpretation in terms of novel forms of nuclear matter 
such as the quark-gluon plasma. 
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